Topics Concerning the Quadrupole-Quadrupole Interaction 
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J^j ■ Abstract 
O 

ON . 

We address some properties of the quadrupole-quadrupole (Q ■ Q) interaction 
in nuclear studies. We first consider how to restore SU(3) symmetry even 
though we use only coordinate and not momentum terms. Using the Hamil- 

3 : tonian H = + f<A ? ) " xEi< 7 Q« • Qti) - I EiC?« • Q(i) with 



= r 2 Y2^, we find that only 2/3 of the single-particle splitting (eod — ei s ) 
comes from the diagonal term of Q ■ Q -the remaining 1 /3 comes from the in- 
teraction of the valence nucleus with the core. On another topic, a previously 
derived relation, using Q-Q, between isovector orbital B{M1) (scissors mode) 
and the difference (B(E2,isoscalar) — B(E2, isovector)) is discussed. It is 
shown that one needs the isovector B(E2) in order that one get the correct 
limit as one goes to nuclei sufficiently far from stability so that one subshell 
(neutron or proton) is closed. 
In this work we address issues pertaining to shell model calculations with the schematic 
quadrupole-quadrupole interaction. Even today, this interaction is of value in casting light 
upon the relationship between shell model and collective model behaviour. There are still 



new things to be learnt about this interaction in nuclei, and we will discuss two examples 
here. 



I. THE SINGLE-PARTICLE SPLITTING (e 0D - e is ) NEEDED TO GET THE 577(3) 

RESULT 

We wish to obtain Elliott's SU(3) results Jl| in a shell model calculation in which only the 
coordinate Q ■ Q interaction is used. We do not wish to use the momentum-dependent terms. 
The latter were introduced by Elliott so that, in combination with the coordinate terms, 
there would be no AiV = 2 admixtures i.e. no admixture from configurations involving 2 
hu excitations. However, we want to see the effects of such admixtures in our shell model 
studies. One classic problem in which AiV = 2 admixtures are important id the E2 effective 
charge, but there are many other problems of interest along these lines. 

The Hamiltonian we consider is therefore 

H = + \™ 2 r\) - x E Q(i) ■ Qti) - § E Q® ■ Q® 

i i<j i 

where Q{i) k ■ Q(j) k = (-l) k \/2k + lr(i) k r(j) k [Y(i) k Y(j) k }° with k = 2. Like Elliott, we 
have not only the two-body Q ■ Q term, but also the i = j single-particle term. 

It is convenient to introduce the following quantity: x = 56 4 x/327r where b is the oscil- 
lator length parameter, such that b 2 = h/muj = 41.46/ hu>. 

To evaluate the single-particle term we use the addition theorem: 

P k (cos9 12 ) = ^i^£y^(l)>V M (2) 

thus 

The single-particle potential is then 

U(r) = • == 4x(Q 4 
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The expectation value of this single-particle term for various single-particle states is given 
in Table I. What single-particle splitting e od — is needed to get Elliott's SU(3) results? 
The best way to answer this is to give the formula for the 577(3) energy in the Is — Od shell 
(in which the momentum terms are included): 

E(\fj) = x [-4(A 2 + /i 2 + Xfi + 3(A + //)) + 3L(L + 1) 

For a rotational band, the L = 2 — L = splitting is given by the last term and is equal to 
18%. This must also be the e d — €\ s because it is also an L = 2 — L = splitting. But, as seen 
from Table I, the splitting due to the diagonal Q ■ Q interaction is (—63 — (—75))% = 12%. 
Where does the remaining 6% come from? 

The answer is that the missing part comes from the interaction of the particle with the 
core. For Q ■ Q, the only contribution is the exchange term of the Od particle with the Os 
core. 

Thus, to get the Elliott S77(3) results in the Is — Od shell, we must not only include his 
diagonal term but also include the particle-core interaction i.e. take the shell model as an 
A particle problem rather than an (A — 16) particle problem. 

The same thing happens in the 0/ — lp shell. The single-particle splitting required to 
get the SU(3) result is e / — e\ v = 3(3 x 4 — 1 x 2)x = 30%. As seen from Table I, we only 
get 2/3 of this (20%) from the diagonal Q ■ Q term. The remaining 10% comes from the 
interaction of the valence nucleons with the core (actually only the Op shell in the core will 
contribute) . 

II. CLARIFICATION OF A RELATION BETWEEN THE ISOVECTOR ORBITAL 
MAGNETIC DIPOLE TRANSITION RATE (I.E. SCISSORS MODE EXCIATION 
RATE) AND THE ELECTRIC QUADRUPOLE TRANSITION. 

As a second example, we will attempt to clarify a relationship between orbital magnetic 
dipole transition rates (i.e. scissors mode excitation rates) and electric quadrupole transi- 
tions. Using the interaction — xQ • Q, Zheng and Zamick H obtained a sum rule relating 



these two quantities. The isovector orbital magnetic dipole operator is (L n — L u )/2 (the 
isoscalar one is half the total orbital angular momentum L/2 = (L n + L u )/2). In detail, the 



sum rule reads 
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where B(M1) is the value for the isovector orbital Ml operator (g^ = 0.5 g\ v = —0.5 
9 an = g su = 0) and the operator for the E2 transitions is J2 pr otons e P ^ 2 ^2 + Uneutrons ^nr 2 Y 2 
with e p — 1, e n — 1 for the isoscalar transition (IS), and e v — 1, e n — —1 for the isovector 
transition (IV). The above result holds also if we add a pairing interaction between like 
particles i.e. between two neutrons and between two protons. 

The above work was motivated by the realization from many sources that there should 
be a relation between the scissors mode excitation rate and nuclear collectivity. Indeed, the 
initial picture by Palumbo and Loludice || was of an excitation in a deformed nucleus in 
which the symmetry axis of the neutrons vibrated against that of the protons. In a 1990 
contribution by the Darmstadt group |5|, it was noted that the Sm isotopes, which undergo 
large changes in deformation, the B(Ml) SC i SSOrs was proportional to B(E2, 0\ — > 2i). The 
B(E2) in turn is proportional to the square of the nuclear deformation S 2 . 

The above energy-weighted sum rule of Zheng and Zamick was an attempt to obtain 
such a relationship microscopically using fermions rather than interacting bosons. To a 
large extent they succeeded, but there are some differences. Rather than the proportionality 
factor B(E2,0i — > 2i), there is the difference of the isoscalar and isovector B(E2). Now 
one generally expects the isoscalar E2 state to be most collective and much larger than the 
isovector B(E2). If the latter is negligible, then indeed one basically has the same relation 
between scissors mode excitations and nuclear collectivity, as empirically observed in the 
Sm isotopes. 

However, derivation of the above energy-weighted sum rule is quite general, and should 
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therefore hold (in the mathematical sense) in all regions -not just where the deformation is 
strong. To best illustrate the need for the isovector B(E2), consider a nucleus with a closed 
shell of neutrons or protons. In such a nucleus, and neglecting ground-state correlations, the 
scissors mode excitation rate will vanish -one needs both open shell neutrons and protons 
to get a finite scissors mode excitation rate. On the other hand, the B(E2, Oi — > 2i) can be 
quite large. However, if we have say an open shell of protons and a closed shell of neutrons, 
the B(E2,0i — > 2i) can be quite substantial. Many vibrational nuclei are of such an ilk, 
and they have large B(E2)'s from ground e.g. 20 W.u. 

However, in the above circumstances, the neutrons will not contribute to the B(E2) 
even if we give them an effective charge. But if only the protons contribute, it is clear that 
B(E2, isovector) = B(E2,isoscalar). 

As an example, let us consider the even-even Be isotopes 6 Be, 8 Be, 10 Be and 12 Be. In so 
doing, we go far away from the valley of stability, but this is in tune with modern interests 
in radioactive beams. 

Fayache, Sharma and Zamick [[J have previously considered s Be and 10 Be. The point 
was made that these two nuclei had about the same calculated B(E2,0 1 — > 2!), but the 
isovector orbital 5(Ml)'s were significantly smaller in 10 Be than in 8 Be. This was against 
the systematic that B{Ml) or utai is proportional merely to B(E2). In detail, the calculated 
B(M1,0 1 -> U) was 2/tt/4 for 8 Be, and in w Be it was 9/32tt/4 (T = 1 -> T = 1) and 
15/32tt/4 (T = 1 -»• T = 2) . Thus the ratio of isovector orbital B(Ml)'s 10 Be/ 8 Be = 3/8. 

We now extend the calculations to include 6 Be and l2 Be. These are singly closed nuclei. 
We see in Table II how everything hangs together. We can explain the reduction in B(M1) 
in 10 Be relative to 8 Be by the fact that the isovector B(E2) in 10 Be is much larger than in 
8 Be. Note that the isoscalar B(E2)'s are almost the same in these two nuclei. The summed 
B(M1) in 8 Be is 2/tt, but in 10 Be it is only 3/8 of that. 

In 6 Be and l2 Be, the E2 transition is from two protons with L = S = to two 
protons with L = 2 S = 0. Note that, surprisingly, the coefficients in front of the effective 
charge factors is larger for singly-magic 6 Be than it is for the open shell nucleus 8 Be. The 
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factors are respectively 12.56 4 /7r and 8.756 4 /7r. However, the charge factor for G Be ( 12 Be) 
is e 2 , whereas for s Be it is (e p + e n ) 2 . The latter gives a factor of four enhancement for the 
isoscalar B(E2) in 8 Be. 

Again we see from Table II that the isoscalar and isovector B(E2)'s are necessarily the 
same and, when this is fed into the sum rule of Zheng and Zamick one gets the consistent 
result that B(Ml) or utai is zero. 

This work was supported by a Department of Energy grant DE-FG02-95ER 40940. One 
of us (L.Z.) thanks E. Vogt, B. Jennings, and H. Fearing for their help and hospitality at 
TRIUMF. 
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TABLES 

TABLE I. The Expectation Value of U/x = —^{r/b) A for several single-particle states 

state (U/x) 

Os -15 

Op -35 

Od -63 

Is -75 

0/ -99 

lp -119 



TABLE II. The Values of B{Ml) orhita i and B(E2) isoscalar and B(E2) isosvector for Be isotopes. 

Nucleus B(Ml) orbltal B(E2) isoscalaT (e 2 /m 4 ) a B{E2) isosvector {e 2 fm A f 

6 Be 19.92 b 19.92 b 

8 Be 2/vr=0.637 82.66 c 7.371 

10 Be T = 1 — ► T = 1 9/32^=0.0895 70.78 d 32.35 e 

T = 1 -► T = 2 15/32tt=0.149 3.322 

12 5e 19.92 19.92 



a The value b = 1.650 fm was used for all nuclei above. 

b The analytic expression in 6 Be is B{E2) = £6 4 e 2 . 

c The analytic expression in 8 Be is B(E2) = ff & 4 (e p + e„) 2 . 

d B(E2) isosca [ ar = to the 2^ state, and is equal to 68.24 e 2 fm A to the 2% state. 

c B(E2)i sovector = 31.19e 2 /m 4 to the 2^ state, and is equal to zero for the 2^ state. 
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